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The most precise comparison between theory and experiment for the B Xst^i~ rate is in the 
low region, but the hadronic uncertainties associated with an experimentally required cut on mx 
potentially spoil the search for new physics in these decays. We show that a 10-30% reduction of 
dr(B ^ Xst*'l~)fdq^ due to the mx cut can be accurately computed using the B Xs"f shape 
function. The effect is universal for all short distance contributions in the limit m^ <C m%, and 
this universality is spoiled neither by realistic values of the mx cut nor by as corrections. Both the 
differential decay rate and forward-backward asymmetry with an mx cut are computed. 


I. INTRODUCTION 


In the standard model (SM) the flavor-changing neu¬ 
tral current process B does not occur at tree 

level and is thus a sensitive probe of new physics. Pre¬ 
dicting its rate involves integrating out the W, Z, and t 
at a scale of order by matching on to the Hamilto¬ 
nian ill, 01 


Hw = -^VtbV* 


6 10 


i=7 


, ( 1 ) 


evolving to fi = frih, and computing matrix elements of 
Hw ■ Here Oi — Oe are four-quark operators and 


been studied only in the Fermi-motion model |]rT|. [The 
high region is unaffected by the mx cut, but the 
rate is lower, and calculating it involves an expansion 
in Aqcb/( mb - y^)-] 

In this letter we compute the B Xsi^£~ rate with 
an mx cut in the low region in a model-independent 
framework. For = ©(AgcoWh), the local OPE 

used in all earlier analyses breaks down and must be re¬ 
placed by an OPE involving b quark distribution func¬ 
tions (shape functions), as explained below. We will com¬ 
pute 

f<l2 d^r ■ 


O7 = mb sa^^.eF^^'^PRb, 

Os = rrib PrO, 

O 9 = e 2 (s 7 ^Pi 6 )(£V^), 

Oio = e^(s7^PL6){£7'"75£), (2) 

where Pi,/} = (l=F 75 )/ 2 . Measurements of 07 , 8 , 9.10 probe 
flavor-changing neutral currents and test the SM. This 
can be done with the dilepton invariant mass spectrum, 
dV/dq^, with q^ = {p^+ It is calculable in an 

operator product expansion (OPE), and the nonpertur- 
bative corrections are 0 (Aq(;,q/to^) [§, jj]. The matching 
and anomalous dimensions for Ci are known at next-to- 
next-to-leading log (NNLL) order, as are the perturbative 
QCD corrections to the matrix elements of Oi |||, ||, 
(except small Oa-e terms). 

A complication in P ^ Xsi^(.~ compared with B 
XsU is that the long distance contributions, B J/ipXg 
and tp'Xs followed by J/ilt, ijj' are two orders 

of magnitude above the short distance prediction, a fact 
which is not well understood. Therefore, either theory 
and data are both interpolated, or the short distance 
calculation is compared with the data for q"^ < 

or q'^ > . The low q^ region, q^ < 6 GeV^, al¬ 

lows the most precise comparison with the SM, but re¬ 
quires a cut on the invariant mass of the hadronic fi¬ 
nal state, mx < m™*. In the latest Belle analysis 
= 2GeV §, while Babar uses = 1.8 GeV §. 
This cut is to remove backgrounds and will likely be re¬ 
quired for quite some time [^. So far, its effect has 


and study the ratios 


Vij 91; 92 ) 


pcut 



( 4 ) 


For convenience we define normalization factors 
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with kinematic dependence G 99 = Gqo = (2q^ -I- ml), 
G 77 = Am'^p{l + 2m1/q^), and G 79 = 12mBmb. Here and 
below, mb is a short distance mass, such as ml^ |Q. In 
Eqs. (§-^, ij = {77, 99, 00, 79} label contributions of 
time-ordered products of operators, r{Oj,Oi}. The to¬ 
tal decay rate with cuts is the sum of these contributions. 


pcut 
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mx 

qf<q^<qi 


Err- 
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We will also study 77 / = 77 / ( 71 }}'^“*, gj, g^), which differ 
from r/ij by the replacement of mx by pj} = Ex — \px \ '■ 
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In Eqs. (^-|^) the short distance coefficients 07 , 9,9 track 
the G 7 , 9 ,io dependence in Eq. (||) that one would like to 
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measure. Here c^ = cg = and cg = 

Cio can be obtained from local OPE calculations Q at 
each order, as discussed in Ref. Q. 

The ? 7 ij’s contain the effects of the mx cut, and are 
defined with a normalization that makes them less sensi¬ 
tive to the choice of ij. At leading order in Aqcd/w;, and 
as, rjij give the fraction of events with mx < and 

r]ij = 1 for = ms- This interpretation is altered at 
subleading order by as corrections, but knowing rjij at a 
given order in perturbation theory is still sufficient to de¬ 
termine r™* and thus the total rate with cuts in Eq. 
at this order. In principle, depend in a nontrivial 
way on ij (and qf and due to different dependence 
on kinematic variables, Og corrections, etc. At leading 
order in Agco/'m-bj we demonstrate that rjij are actually 
independent of the choice of ij, a property which we call 
“universality”. We first show this formally in section II 
at leading order in p~^/ms ^ 1 for Then in 

section III we demonstrate it numerically for the exper¬ 
imentally relevant 77 ( 771 ™*), including the as corrections 
and phase space effects. 

We also compute the decay rate for current experimen¬ 
tal cuts. We find that the rate is sensitive to the choice 
of the cut, and that the cut causes a reduction in the 
rate by 10-30%. Since the same shape function occurs in 
B Xs£~^£~, Xu£p, and Xsj, the tti™* or Px'^'^^ depen¬ 
dence in one can be accurately determined from the oth¬ 
ers, and the magnitude of the reduction can be computed 
quite accurately. Alternatively, instead of using the theo¬ 
retical computation of the 7773 )** dependence, universality 
can be exploited to remove the main uncertainties, by 
normalizing the B Xs£~^£~ rate to B —!■ Xu£i'- 


II. mx CUT EFFECTS AT LEADING ORDER 


For simplicity, consider the kinematics in the B me¬ 
son’s rest frame. Since q = Pb — Px, 

2mBEx = rn% + m\ - q^. ( 8 ) 

If m\ Wb and q^ is not near 777 ^, then Ex = 0{mB)- 
Since E\ ^ 777 ^, px is near the light-cone, with = 
Ex - \px\ = O(Aqcd) and p-^ = Ex + \px\ = 0{mB)- 
Of the variables symmetric in p^+ and pt- (p^, Ex, q^, 
m\), only two are independent, and we work with q^ and 
pj3 or mx- The phase space cuts are shown in Fig. |^. 

For the pj). <C p^ region, factorization of the form 
dr = E[J®f^^'> has been proven for semileptonic and 
radiative B decays | p^ , where H contains perturbative 
physics at Pb ^ 'mb, J at pt ~ y^AqcDWb, and 
is a universal nonperturbative shape function llffl. This 
factorization also applies for B Xs£^£~ with the same 
/(O) 

, as long as q^ is not parametrically small 0- 
In the q^ < 6GeV^ region, \C^^^{q^,po = 4.8GeV)| = 
4.52 to better than 1%, and can be taken to be constant. 



FIG. 1: Phase space cuts. A substantial part of the rate for 
qi < < qi falls in the rectangle bounded by < pj'’”*. 


We neglect as corrections in this section and find 

dr ^ Ao)/ -t') Zo_ [(777 b -pj3)^ 
dpjdq2 m% - p+y 

X {(|Cg“-|2 + C?o)[2g2 + (mB-pi)^] 


■4777b \C^ 


mix 1 2 


1 + 


2 ( 777 b - Px) 


,+ f2. 


+ 12r77BRe[Cr"C9““*](777B -p^)}, 


(9) 


where has support in w G [0, 00 ). As a function 

of p^, the kinematic terms in Eq. (^ vary only on a 
scale 777 b, while varies on a scale Aqcd- Writ¬ 

ing 777b = 777f, -I- A and expanding in (pj^ — A) /tub decou¬ 
ple the Px and q^ dependences in Eq. (ffl, and give ex¬ 
actly the local OPE prefactors, {ml — Gij{q^), used 
in Eq. (|). For 77^(pi™‘, 9i, g|), the pj^ integration is 
over a rectangle in Fig. whose boundaries do not cou¬ 
ple pj3 and q^. Thus, with the above expansion, we find 
p'ij = rj’, where 


= j dpi . 


( 10 ) 


independent of ij and q^, q^. While the mx cut retains 
more events than the p^ cut, the latter may give theoret¬ 
ically cleaner constraints on short distance physics when 
statistical errors become small. 

The effect of the mx cut is dependent, because the 
upper limit of the p^ integration is dependent, as 
shown in Fig. 1. When we include the full p^ depen¬ 
dence in Eq. (3), the universality of rjijim'^x^,q\,qD is 
maintained to better than 3% for 1 GeV^ < g? < 2 GeV^, 
5GeV^ < gi < 7GeV^ and 7775 ^“ > 1.7 GeV, be¬ 
cause the region where the p(^ and g^ integration lim¬ 
its are coupled has a small effect on the ij dependence. 
This is exhibited in Fig. ||, where the solid curves show 
Pij{mPx^, 1 GeV^, 6 GeV^) with the shape function set to 
model 1 of Q, with ml^ = 4.68GeV and Ai from [^ . 
(Taking qf = 1 GeV^ instead of increases the sensi¬ 
tivity to Ggpo, but one may be concerned by local dual¬ 
ity / resonances near g^ = 1 GeV^. To estimate this un¬ 
certainty, assume the (j) is just below the cut and B{B 
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FIG. 2: ^7iJ(m5J'^lGeV^6GeV^) as functions of . The 
dashed curves show the local OPE result, the solid curves 
include the leading shape function effects. The uppermost, 
middle, and lowest curves are r;oo, 99 , f? 79 , and 7777 , respectively. 


Xscf) « 10 X B{B Then B Xs 4 > X.t^l- 

is ~2% of the Xst^l~ rate.) 

The local OPE results for f?!) obtained 

by replacing f^^\px) by 6 {A—px) in Eq. (^. Performing 
the Px integral sets [ms — p\) = and implies m\ > 
A(mB — 9^/tnt,). This makes the lower limit on equal 
max{g^, mb[mB — (w™*)^/A||, and so the rjij's depend 
on the shape of dPij. In Fig. ||the local OPE results are 
shown by dashed lines, and clearly 7777 ^ 7799. However, 
the local OPE is not applicable for ^ Aqcd- 

The universality of 77^ found here could be broken by 
Us corrections in the H or J functions, or by renormaliza¬ 
tion group evolution, since these effects couple Px and 
and have been neglected so far. We consider these next. 


III. CALCULATION AND RESULTS AT 0{as) 

A complication in calculating B Xsi^i~ compared 
with B Xuf.v is that, in the evolution of the effective 

I 


Hamiltonian down to mb, Cq{p) receives a ln(77i^/m^) 
enhanced contribution from the mixing of 02 - Thus, for¬ 
mally, Cg ^ 0 {l/as), and conventionally one expands 
the amplitude in as, treating as In(777^/777^) = 0 {\) [^ . 
In the local OPE this is reasonable, since the nonpertur- 
bative corrections are small, and at next-to-leading log 
(NLL) all dominant terms in the rate are included. How¬ 
ever, in the shape function region nonperturbative effects 
are 0(1) and only the rate is calculable, not the ampli¬ 
tude. With the traditional counting, the Cg contribution 
to the rate would be needed to 0(as) before the CIq 
terms could be included. 

This would be a bad way to organize the perturbative 
corrections (numerically \Cg(mb) \ ~ lOiol). It can be cir¬ 
cumvented by using a “split matching” procedure to de¬ 
couple the perturbation series above and below the scale 
mb 0 . This allows us to consider the short distance 
coefficients C™, ( 7 ™“, and Cio as 0 {\) numbers when 
organizing the perturbation theory at m^ and mb Aqcd ■ 

The rate and the forward-backward asymmetry are 


d^r 

dgMpx 

d^ApB 

d<z2dp+ 


^H{q^,p-^)F^°\p+,p ), 
mB 

^ K{q^,p'l.) F^^\p+,p-) , 
rriB 


( 11 ) 


where p~ = mb — q^ j (tub ~ Px)- The hard functions 
FI and K were computed in using soft-collinear ef¬ 
fective theory (SCET) ^o| and split matching. This 
factorizes the dependence on scales above and below mb 
as Tij ^ F[i{iJ,Q)H2{pb)F^^\pb), with separate pg and pb 
independence. Up to the order one is working at, Fli is pg 
independent, the pb dependence in Hg and F^°^ cancels, 
and is Pi independent. The shape function model is 
specified at pA- The convolution of jet and shape func¬ 
tions at NLL including as corrections is 


F^°\Px>P ) 




Uh{p ,Pi,Pb) |^21n^^^f- 31 n^^|--h 7 - 




+ jg Y ^^p^ ~ (Pi(l - {ph Pi)) 

gVs(Mi.MA) / uj /■! , , 

— ) / dt/(°)[w(l-t^/’'),/7A], 


r(i + ?7) \phj 


( 12 ) 


where 77 = ( 16 / 25 ) ln[as(/7A)/as(/ii)], Uh was computed in Ref. |^, the one-loop jet function in Ref. H^, ^ , and the 
shape function evolution up to pi in Refs. ||^, (for earlier calculations, see Refs. 0il)- The hard coefficients 
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H and K for B ^ Xgt^l are 
.2 


H{q ,Px) = 


K{q\p-^x) = - 


(l-P+)3 

+ A\Cf^{po)\^[nUs,Pi,) + 

3g2[(l -p+)2 _ q2]2 


2(l-p+)2 


(l-P+)3 

where s = q^jml, (p = q^jm^, p~^ = p'^/ttib, and 


■ ri|,(s, pb) + 12Re[C“‘’'(^o)C'““(s, ^ 0 )*] (1 - Px)^Eis, pt) 

y 

^2^(s,/if,)Re{c*o[c9”“(s,Mo)f^A(s,Mfa) + pb)]}, (13) 


a.c 


— 1 H-( 5 , fib) j — 1 H-( 5 , fib) + 


TT 

TT 


TT L 

a, 

TT 


(l-p+)2-^2 

2 (l-p +)2 


^bis)+ujY{s) 


ric = f +-^ u;^{s, Pb), = 1 + ^ [wj(s, ^h) — a;J(s)] , + llsllc) ■ (14) 


Here as = as{pb) and are defined in Ref. Il4| . 

For each shape function model, the deviations of the 
Pij's from being universal, with all NLL corrections, are 
still below 3%. Thus, the picture of universality in Fig. 
remains valid at NLL order. For this reason we can ex¬ 
plore the overall shift by just studying 7700 . 

In Fig. ^ we plot 7700 (^ 5 ^*, 1 GeV^, 6 GeV^), includ¬ 
ing the as corrections. We use ten different models for 
Our base model has five parameters, three of which 
are chosen to obey the known constraints on its mo¬ 
ments 1 ^, converted to the IS* mass scheme used here. 
For each of five different choices of the remaining two pa¬ 
rameters, we choose two values of the scale, pA, where 
the model is specified. The choice of these ten models 
is guided by making them consistent with the B Xg^ 
data. The ten orange, green and purple (medium, light, 
dark) curves correspond to ml^ = 4.68 GeV, 4.63 GeV, 
and 4.73 GeV, respectively, with the central values po = 
Pb = 4.8 GeV and pi = 2.5 GeV. The curves with slightly 



FIG. 3: 7700 ( 777 ^*, 1 GeV3, 6 GeV^) as a function of 7773 ^**. The 
orange, green and purple (medium, light, dark) curves show 
mP = 4.68 GeV, 4.63 GeV, and 4.73 GeV, respectively. 


lower [higher] values of poo at large 777 ™* correspond to 
Pa = 1.5 GeV [2 GeV]. The spread in the curves gives our 
determination of the uncertainty from the choice of shape 
function model and from mb- For 777 ™* = 2 GeV, varying 
Pb in the range 3.5 GeV < pb < 7.5 GeV changes 7700 by 
±6%. We find a ±5% variation for 2 GeV < pi < 3 GeV. 

Using the c^’s at NLL, for 1 GeV^ < q^ < 6GeV^ we 
obtain cut branching ratios 

(1.20 ± 0.15) X 10-3 [ 7775 )“ = 1.8 GeV], 

(1.48 ± 0.14) X 10-3 [ 7775 ^* = 2.0 GeV], 

(15) 

where uncertainties are included from nib, Pb, Pi, and 
/(3). Changing po to 3.5 GeV (10 GeV) changes both of 
these rates by —2% (-1-7%), and this uncertainty will be 
reduced by including NNLL corrections i 10- 

The largest source of universality breaking in the 77 ^ ’s 
and one of the largest uncertainties in the cut rate is 
due to subleading shape functions, which affect the rate 
by - 5% for = 2 GeV and by -10% for 7775 )“ = 
1.8 GeV ||. 

If the 777™* dependence were not universal, it would 
modify the zero of the forward-backward asymmetry, 
^fb(9o) = 0- Foi' w™* = 2 GeV we find at NLL 
A< 7 g « —0.04GeV^, much below the higher order un¬ 
certainties i 0- Howeve^ we obtain gg = 2.8 GeV^, 
lower than earlier results [g. The reason is that in the 
SCET calculation of Hfb, using K in Eq. (p^, the pole 
mass 777 ])°^® never occurs, only — p^ and rUb (at this 
order, C™ = {rnb/mB)C^^ @)- Thus, schematically, 
gg ~ 2777b[mf,(/7o)C'f®(77o)]/Re[C'|®(go)], and there is no 
reason to expand nib in terms of 777 ^°*®. 

In the above analysis, the nonperturbative shape func¬ 
tion /(3) was extracted from moments and the B Xgj 
energy spectrum, and this was used as input in deter¬ 
mining our B Xst^t~ results. The overall 10% theo¬ 
retical uncertainty in this approach could be reduced by 
raising the 777 ™*. An alternative approach would be to 
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keep TO™* < rnu and measure 


r™*(B ^ Xsi+i-) 

r™*(s ^ XJp) 


(16) 


with the same cuts used in the numerator and denomina¬ 
tor. The dependence of the semileptonic rate on to™* is 
identical to that of Tqq*. A measurement of R bypasses 
the need for a shape function model, because we found 
that the mx-cut effects are universal to a very good ap¬ 
proximation and therefore cancel between the numerator 
and denominator of R. 

In conclusion, we pointed out that the observed B —> 
Xsi~^i~ rate in the low region is sensitive to the exper¬ 
imental upper cut on mx ■ The reduction in the rate due 
to this cut is determined by the universal B shape func¬ 
tion. In the region of the experimental measurements 
an OPE exists only for the decay rate and not for the 
amplitude, a fact that necessitates a reorganization of 


the usual perturbation expansion. Since one can use the 
shape function measured in other processes, the sensitiv¬ 
ity to new physics is not reduced. We found that the 
77 ’s for the different operators’ contributions are univer¬ 
sal to a good approximation. These results also apply for 
B —> Xd£~^i~, which may be studied at a higher lumi¬ 
nosity B factory. Subleading Aqcd / mb as well as NNLL 
corrections to the rate and the forward-backward asym¬ 
metry will be studied in a separate publication [ p4| . 
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